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ABSTRACT 
Measurements  of  the  ultrasonic  attenuation  of  10  and  30  Mcs/sec 
longitudinal  waves  by  pulsed-echo  techniques  were  made  on  a  99.999% 
pure  single  crystal  of  superconducting  zinc  in  the  0001   direction 

as  a  function  of  temperature  from  4.2  K  to  0.320  K  using  an  open-ended 

3 
type  He   cryostat.   The  attenuation  was  found  to  be  frequency  dependent 

and  decreased  less  sharply  near  the  superconducting  transition  tem- 
perature, T  ,  than  predicted  by  the  Bardeen-Cooper-Schrief f er  (BCS) 
theory.  Attenuation  due  to  electron-phonon  interactions  only  was 
found  by  subtracting  from  experimental  points  the  value  of  residual 
attenuation  gotten  by  extrapolation  of  the  data  to  T  =  0  K.  Using 
the  BCS  theory  the  zero  degree  superconducting  energy  gap  was  found 
to  be  2  £(0)  =   (3.36  +  0„13)kT  with  T  =  0.817°K. 
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INTRODUCTION 


The  rapid  decrease  in  the  attenuation  of  ultrasonic  waves  in 

superconductors  below  the  superconducting  transition  temperature,  T   , 

12-  ° 

was  first  observed  by  Bommel  and  Mackinnon.  '     Bommel  investigated 

the  ultrasonic  attenuation  of  shear  and  longitudinal  waves  in  a  cylin- 
drical crystal  of  lead  in  the  frequency  range  from  9-27  Mcs/sec  from 
approximately  18  K  to  1.5  K.   He  found  a  marked  decrease  in  the  atten- 
uation in  the  superconducting  state  which  was  frequency  dependent.   No 
theoretical  explanation  was  made,  though  he  suggested  that  the  results 
might  lend  support  to  the  idea  that  superconductivity  is  due  mainly  to 
interactions  between  electrons  and  lattice  vibrations.   Mackinnon  in- 
vestigated the  absorption  of  10  Mcs/sec  sound  pulses  in  a  superconducting 
polycrystalline  tin  rod  and  in  a  polycrystalline  magnesium  rod  between 
1.5  K  and  4.2  K.   In  tin  he  found  the  attenuation  decreased  by  a  factor 
of  about  50  and  in  magnesium  that  the  attenuation  decreased  steadily 
with  decreasing  temperature. 

This  decrease  in  the  electron  contribution  to  the  ultrasonic  atten- 
uation was  first  satisfactorily  explained  in  1957  by  the  theory  of 

3 
Bardeen,  Cooper,  and  Schrieffer  (BCS  theory).    The  BCS  theory  predicts 

that  the  ratio  of  the  electron  contribution  to  the  attenuation  in  the 

superconducting  state,  a      ,  to  that  in  the  normal  state,  ot      ,  is 

s  n 


(X 

s 

2 

n 

exp 

r«<r)i 

+ 1 

where   2e(T)  is  the  temperature  dependent  energy  gap.   The  theory  further 
predicts  that  as  T  approaches  zero,  the  energy  gap  approaches  3.50  kT  . 

4 
Morse  and  Bohm  were  the  first  to  apply  this  theory  to  the  analysis 

of  ultrasonic  attenuation  data  on  tin  and  indium,  shortly  before  the 

publication  of  the  BCS  theory,  basing  their  analysis  on  a  conversation 


with  Professor  Bardeen .   Their  analysis  pointed  out  the  remarkable 
agreement  between  the  BCS  theory  and  their  experimental  data.   Since 
then  the  BCS  theory  has  been  applied  to  many  superconducing  metals.  ' 

These  investigations  have  all  shown  close  agreement  between  the 
theory  and  experiment,  but  have  also  indicated  the  need  for  refinement 
of  some  of  the  simplifying  assumptions  of  the  BCS  theory.   The  theory 
shows  that  the  basic  interaction  causing  superconductivity  is  that  of 
the  exchange  of  a  virtual  phonon  between  two  electrons.   An  electron 
moving  through  the  lattice  creates  a  virtual  phonon  which  is  absorbed 
by  a  nearby  electron.   If  the  distortion  of  the  lattice  is  such  as  to 
surround  the  electron  by  a  positive  charge  greater  than  the  electronic 
one,  the  second  electron  sees  an  attractive  potential  greater  than  the 
repulsive  Coulomb  one  and  is  attracted.   These  electrons  form  what  is 
known  as  Cooper  pairs.   The  theory  then  shows  that  the  lowest  energy  of 
the  system  is  achieved  by  the  maximum  pairing  of  electrons  of  opposite 
spin  such  that  their  momenta  add  to  zero.   The  difference  in  the  normal 
ground  state  energy  and  the  superconducting  ground  state  energy  is  due 
entirely  to  an  interaction  potential,  less  than  zero,  arising  from  the 
correlation  between  the  Cooper  pairs.   The  energy  of  this  interaction 
leading  to  the  scattering  of  one  pair  to  another  is  the  difference 
between  the  energy  of  interaction  by  means  of  a  phonon  and  the  energy 
of  Coulomb  interaction. 

The  existence  of  an  energy  gap  then  follows  from  the  small,  but 

finite,  amount  of  energy  required  to  break  up  a  Cooper  pair.   But  more 

important,  having  a  single  unpaired  electron  removes  a  large  number 

of  pairs  which  could  have  occupied  the  state  of  a  Cooper  pair.   This 

gap  will  be  temperature  dependent  since  increasing  numbers  of  pairs 

will  be  broken  up  by  thermal  excitation  as  the  temperature  rises  from 

0°K  to  T  . 
c 

Morse,  Olsen,  and  Gavenda  pointed  out  the  need  to  consider  the 
anisotropy  of  the  energy  gap  and  showed  that  in  tin  there  is  considerable 


anisotropy  (3.2  kT  -  4.32  kT  ) .   In  order  to  conserve  both  momentum  and 
r-/        c  c 

energy  in  a  phonon-electron  interaction,  the  component  of  the  electronic 
velocity  parallel  to  the  direction  of  sound  must  equal  the  phonon  veloc- 
ity; i.e.,  the  sound  velocity.   Since  the  velocity  of  electrons  at  the 
Fermi  surface  is  in  general  much  greater  than  that  of  sound,  only  those 
electrons  moving  at  essentially  right  angles  to  the  direction  of  sound 
propagation  can  interact  with  the  phonons .   On  the  Fermi  surface  these 
electrons  lie  in  a  narrow  band  perpendicular  to  the  sound  velocity. 
Since  the  ultrasonic  attenuation  in  a  particular  direction  is  an  average 
over  this  narrow  band,  the  energy  gap  measured  will  be  an  average  for 
the  electrons  lying  within  the  band.   Anisotropies  in  the  Fermi  surface 
will  then  lead  to  anisotropies  in  the  energy  gap.   The  BCS  theory  makes 
the  simplifying  assumption  that  the  interaction  energy  is  isotropic  and 
constant  near  the  Fermi  surface .   It  is  this  assumption  that  leads  to 
the  BCS  prediction  that  the  zero  degree  energy  gap  is  the  same  for  all 
superconductors  and  is  equal  to  3.5  kT  . 

The  relation  for  the  ratio  of  ultrasonic  attenuation  in  the  super- 
conducting state  to  that  in  the  normal  state  follows  from  the  coherence 
of  paired  electrons.   In  the  normal  state  scattering  of  an  electron  to 
another  state  is  independent  of  any  other  transitions.   However,  the 
pairing  of  electrons  in  the  superconducting  state  makes  these  transitions 
interdependent  and  the  theory  shows  that  the  scattering  by  pairs  is 
destructive.   Thus,  the  absorption  of  phonons,  causing  ultrasonic  atten- 
uation, decreases  sharply  below  T  .   With  the  assumption  of  low  energy 

phonons;  i.e.,  hu>  «   2e(0)  ,  the  ratio  at   /at        follows. 

s   n 

The  present  work  on  zinc,  which  has  a  hexagonal  crystal  structure, 

was  undertaken  to  determine  the  energy  gap  in  the  direction  of  the 

c-axis  and  to  compare  the  results  with  measurements  made  in  the  base 

plane.   Ultrasonic  frequencies  on  the  order  of  10  cps  were  used, 

well  satisfying  the  condition  for  the  given  expression  of     at   /a 

s  n 
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EXPERIMENTAL  PROCEDURES 

A  cylindrical,  single  crystal  of  zinc,  approximately  1\"    long  and 
3/4"  in  diameter  with  the  c-axis,  [0001]  direction,  parallel  to  the 
cylinder  axis  was  purchased  from  Semi-Elements,  Inc .   The  purity  of  the 
sample  was  stated  to  be  in  excess  of  99.999%.   From  this  ingot  the  experi- 
mental sample  was  prepared  to  the  approximate  desired  size,  limited  by 
the  cryostat,  by  cleaving  both  ends  and  etching.   The  ingot  was  held  in 
a  styrofoam  form,  immersed  in  liquid  nitrogen  and  cleaved  by  laying  a 
new,  sharp,  razor  blade  at  the  desired  point  and  rapping  it  with  a 
mallet.   To  reduce  the  diameter,  the  ends  were  coated  with  brushing 
lacquer  and  the  ingot  etched  in  dilute  HC1  acid. 

Standard  ultrasonic  pulse-echo  techniques  were  used  and  will  be 
discussed  only  in  terms  of  the  experiment.   The  success  of  this  method 
depends  strongly  on  having  the  end  faces  of  the  sample  flat,  parallel, 
and  perpendicular  to  the  direction  in  the  crystal  one  wishes  to  investi- 
gate.  To  achieve  this,  the  zinc  sample  was  cast  with  a  cerro-alloy  in 
a  bored  out,  cylindrical  jig  whose  faces  were  machined  flat,  parallel 
and  perpendicular  to  the  crystal  axis  to  a  tolerance  of  0.003".   The 
ends  were  first  lapped  on  metallographic  paper  and  then  polished  on  a 
diamond  dust  wheel. 

The  final  sample  was  (3.1602  +  0.0002)  cm  long  and  approximately 
5/8"  in  diameter  at  room  temperature.   A  Laue  back  reflection  pattern 
of  the  zinc  crystal  showed  concentric  circles,  indicating  a  thin  poly- 
crystalline  layer  on  each  end.   However,  the  x-ray  spots  showing  through 
the  circles  indicated  that  the  faces  were  perpendicular  to  the  c-axis. 

A  %",  X-cut,  10  megacycle  quartz  transducer  was  next  bonded  to  the 
zinc  sample.   The  success  of  this  bond  is  a  critical  step  in  the  experi- 
mental procedure.   After  some  experimenting,  the  following  procedure 
proved  to  be  reliable.  Both  the  quartz  and  sample  faces  are  cleaned  first 
with  acetone  to  remove  greasy  substances,  second  with  alcohol  to  remove 

_4- 


water,  and  finally  with  acetone  to  ensure  the  surfaces  are  dry.   The 
best  bonds  were  obtained  on  days  when  the  humidity  was  low.   Next  a 
very  small  amount  of  analytical  grade  glycerin  was  applied  to  the  zinc 
and  the  transducer  centered  and  squeezed  to  make  the  glycerin  flow  and 
also  to  remove,  air  bubbles.   Glycerin  proved  to  be  the  most  satisfactory 
bonding  material,  though  its  use  required  the  making  of  a  new  bond  for 
each  run.   The  critical  temperature  for  the  bond  was  from  40  K  to  50  K. 
If  the  bond  held  on  being  cooled  through  this  range;  the  bond  was  good 

for  the  run.   However,  as  can  be  seen,  one  has  already  started  to  cool 

4 
with  He   and  if  the  run  must  be  stopped  for  a  broken  bond,  the  delay  is 

both  costly  and  time  consuming. 

A  schematic  for  the  ultrasonic  pulsing  equipment  is  shown  in  Fig.  1. 
Radio  frequency  pulses  (10  Mcs/sec  and  30  Mcs/sec)  of  approximately  two 
microseconds'  duration  were  generated  by  an  Arenberg  Pulsed  Oscillator 
triggered  by  a  Time  Mark  Generator  with  a  repetition  rate  of  1000  cps 
or  100  cps.   The  slower  rate  had  to  be  used  below  0.5  K  to  reduce  heating 
of  the  sample.   The  rf  electrical  pulses  were  simultaneously  fed  to  a 
Type  RM  45A  Tektronix  Oscilloscope,  externally  triggered  by  the  Time 
Mark  Generator,  and  impressed  across  the  piezoelectric  quartz  transducer. 
The  zinc  furnished  the  rf  ground  and  the  oscillator,  tuned  to  the  funda- 
mental or  odd  harmonic  of  the  transducer,  caused  it  to  vibrate  and  send 
sound  waves  through  the  zinc.   The  sound  waves  reflected  off  the  other 
end  return  to  the  transducer  to  produce  a  voltage  across  the  electrodes 
furnishing  the  initial  pulse  and  also  to  be  reflected  back  through  the 
sample.  With  a  pulse  rate  of  1000  cps,  these  reflected  echoes  had 
adequate  time  to  decay  away  between  pulses.   The  voltage  produced  by  each 
returning  echo  was  amplified  and  displayed  in  proper  time  relation  to 
the  initial  pulse  on  the  oscilloscope.   The  change  in  height  of  each  echo 
is  a  measure  of  the  attenuation  in  the  sample  due  to  phonon-electron 
interaction,  impurities,  imperfections  and  in  the  transducer  bonding 
system.   Absolute  attenuations  are  gotten  by  taking  pictures  of  the  scope 
traces  with  a  Polaroid  Oscilloscope  Camera,  Model  196A,  Special,  using 


Via.      1 


3000  speed/Type  47A  Polaroid  film  and  plotting  echo  heights  vs  echo 
number  on  semi- log  paper.   Attenuation  effects  other  than  those  due  to 
phonon- e lectron  interactions  are  subtracted  off  as  discussed  in  the 
treatment  of  the  data. 

Data  were  taken  in  the  temperature  range  from  4.2  K  to  approximately 
0.320  K,   The  temperature  range  was  obtained  in  a  He   cryostat  of  the 
open-ended  type,  a  schematic  of  which  is  shown  in  Fig.  2.   It  is  desirable 
to  cool  the  sample  as  fast  as  possible  to  reduce  cooling  time,  but  this 
is  tempered  by  the  necessity  of  not  straining  the  bond. 

The  sample  was  loaded  in  a  copper  sample  holder,  inserted  in  the 

3  3 

He  cryostat  and  the  He   chamber  flushed  with  dry  nitrogen  and  left  at 

3 

one  atmosphere  of  dry  nitrogen.   The  vacuum  jackets  of  the  He   cryostat 

4 
and  the  He   dewar  were  thoroughly  flushed  with  dry  nitrogen  and  left 

4 
soft  at  approximately  300  microns  dry  nitrogen  pressure.   The  He  dewar 

was  also  flushed  and  left  at  one  atmosphere  dry  nitrogen  pressure.   The 

nitrogen  dewar  was  then  filled  with  liquid  nitrogen  and  the  whole  system 

cooled  to  approximately  77  K  which  takes  about  six  hours. 

3       4 
The  He  and  He   jackets  were  then  pumped  out  until  the  pressure  was 

-  fi  "\ 

less  than  5  x  10  mm  Hg .   Keeping  the  pressure  less  than  this  in  the  He 

jacket  was  essential  to  prevent  thermal  conduction  at  the  lowest  temper- 

4 
atures.  While  this  was  being  done,  the  nitrogen  and  He   traps  were 

3 

filled  and  the  He   gas  circulated  for  approximately  one-half  hour  to 

4 
condense  out  all  impurities  except  He  .   This  was  found  to  be  necessary 

as  impurities  would  freeze  out  in  the  cryostat  and  around  the  sample 

3         4 
holder  causing  throttling  when  pumping  full  out  on  the  He  .   The  He 

impurity  caused  no  trouble  as  thermometry  was  done  using  a  carbon  resistor 

3 
and  not  by  using  the  He  vapor  pressure. 

3  4 

When  the  jackets  were  pumped  out  and  the  He   purified,  the  He  dewar 

3 
was  filled  and  about  eight  mm  Hg  of  He   let  into  the  cryostat;  the  re- 
mainder being  returned  to  the  reservoir.   The  system  was  then  cooled  to 
4.2  K,  taking  about  12  hours.   The  process  could  be  accelerated  by 
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3  4 

softening  the  He   jacket  with  He  ,   However,  achieving  a  satisfactory 

vacuum  after  the  system  was  cooled  down  was  very  difficult  since  helium 

at  that  temperature  is  very  difficult  to  pump.   Without  a  good  vacuum 

the  low  temperatures  desired  could  not  be  reached. 

o  4  3 

Once  4.2  K  was  reached,  the  He  was  pumped  on  and  the  He   let 

into  the  cryostat.   A  minimum  temperature  of  1.09  K  can  be  reached 

3 
this  way,  well  below  the  temperature  at  which  He  at  one  atmosphere 

will  condense;  i.e.,  3.195  K.  When  equilibrium  was  reached,  the  cir- 

3 

culating  pump  was  valved  off  to  pump  on  the  He  .   At  500  micons   the 

3 

booster  pump  was  cut  in  and  pumping  full  out  on  the  He   the  lowest 

temperatures  were  reached . 

Temperatures  were  determined  by  measuring  the  resistance  of  a 
%-watt  Speer  Carbon  Resistor,  nominal  resistance  470  ohms,  10%,   grade 
1002,  in  a  Wheatstone  Bridge  Circuit  employing  a  lock-in  amplifier  to 

measure  the  bridge  unbalance.   The  resistor  was  calibrated  against  the 

4  o  4 

vapor  pressure  of  He  down  to  1.1  K  using  the  "1958  He   Scale  of  Temper- 
ature".  Down  to  the  lowest  temperature  reached,  the  resistor  was  cali- 
brated against  the  magnetic  susceptibility  of  cerium  magnesium  double 
nitrate  (CMN) ,  Ce~Mg  (NO  )   »24H  0  using  a  Cryotonics  Mutual  Inductance 
Bridge.    The  susceptibility  of  CMN  is  strongly  anisotropic  with  the 
susceptibility  parallel  to  the  symmetry  axis  essentially  zero.   To  obviate 
the  necessity  of  packing  single  crystals,  powdered  CMN  was  used.   The 
effective  susceptibility  is  then  approximately  657o  of  the  value  perpen- 
dicular to  the  symmetry  axis,  but  the  gain  in  the  packing  fraction  off- 
sets this  decrease. 

Since  the  temperature  dependent  part  of  the  susceptibility  of  CMN 
perpendicular  to  the  symmetry  axis  follows  Curie's  Law 


% 


T 


the  mutual  inductance  bridge  readings  are  correlated  to  the  temperature 
and  thereby  calibrating  the  carbon  resistor. 
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EXPERIMENTAL  RESULTS 

Measurements  of  the  ultrasonic  attenuation  of  10  and  30  Mcs/sec 

longitudinal  waves  were  made  in  the  [0001]  direction  on  a  single  crystal 

of  zinc  whose  ends  were  slightly  polycrystalline .   Data  were  taken  in 

the  temperature  range  4.2  K  to  approximately  0,320  K  and  in  the  earth's 

magnetic  field.   The  superconducting  transition  temperature  was  determined 

to  be  T  =  0.817  K  and  the  zero  degree  superconducting  energy  gap  was 

measured  as  (3.36  +  0.13)kT  . 

—        c 

To  measure  the  attenuation  as  a  function  of  temperature,  pictures 
of  echo  trains  were  taken  at  different  temperatures  and  the  heights  of 
the  echoes  plotted  on  semi- log  paper  to  obtain  the  attenuation  in  units 
of  attenuation  per  echo.   These  were  plotted  as  a  function  of  temperature 
and  when  extrapolated  to  T  =  0  K  gave  an  estimate  of  the  residual  atten- 
uation.  Subtracting  the  residual  attenuation  from  experimental  points 
gave  the  attenuation  due  solely  to  electron-phonon  interactions.   These 
data  were  then  analyzed  using  the  BCS  function  for  reduced  attenuation 
vs  reduced  temperature 


or*  = 


exp  **"•  +1 


The  details  of  all  these  analyses  are  given  in  Appendices  I  -  III. 
In  order  to  illustrate  the  methods,  Figs.  3-14  graphically  display  the 
results  for  one  set  of  data;  those  of  13  November  1964. 

After  the  transition  temperature  was  determined  and  an  estimate  of 
the  residual  attenuation  made,  computer  program  RESATT  was  used  to 
determine  the  residual  attenuation  and  the  zero  degree  energy  gap. 
Figs.  3-5  show  the  point  plots  of  reduced  attenuation  vs  reduced 
temperature  for  three  slightly  different  values  of  residual  attenuation 
subtracted  from  the  experimental  data.   By  inverting  the  BCS  expression 
the  following  function  Y   is  defined: 

Y  =  l/ln(2/a*  -  1)  =  (kT  /e(t))t. 

c 
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Figs.  6-8  depict  the  point  plots  of  the  function  Y  vs  the  reduced 
temperature,   t  ,  corresponding  to  the  same  values  of  residual  atten- 
uation.  All  these  figures  have  the  corresponding  curve  plots  for  the 
BCS  expression  for   2e (0)  =  3.50  kT  . 

The  criterion  used  to  find  the  residual  attenuation  and  the  zero 
degree  gap  was  that  on  the  plot  of  Y  vs   t   the  curve  best  fitting  the 
experimental  points  should  become  a  straight  line  passing  through  the 
origin  for   t   less  then  0.35.   This  follows  from  the  expression  for   Y 
which  becomes  linear  for  2/0f*  much  greater  than  1.00.   The  results  of 
these  analyses  are  given  in  Table  I. 

Once  the  residual  attenuation  was  determined,  computer  program 
VARIGAP  was  used  to  show  how  changing  the  zero  degree  energy  gap  would 
change  the  BCS  curves  for  the  residual  attenuation  and  Y.   Figs,  9-11 
show  the  BCS  curve  for  reduced  attenuation  vs  reduced  temperature  for 
a  given  residual  attenuation  and  with  changes  of  0.10  kT   in  the  energy 
gap  determined  from  RESATT.   Figs.  12-14  give  the  corresponding  results 
for  Y  vs   t .   As  can  be  seen  the  BCS  curves  do  not  shift  much  for  the 
value  of  incremental  change  used. 

It  was  found  that  the  attenuation  was  constant  above  T  and  frequency 
dependent.   Table  I  shows  that  attenuation  increased  with  increasing 
frequency.   The  higher  attenuation  for  30  Mcs/sec  in  the  normal  state 
on  5  January  1965  compared  to  that  of  13  November  1964  reflects  the 
greater  attenuation  losses  from  the  poorer  bond.   The  value  for  the 
zero  degree  gap  was  gotten  from  computations  based  on  semi- log  plots. 
It  is  felt  that  this  method  is  better  than  that  of  comparing  echo 
heights  since  the  latter  method  involves  propagating  through  all  the 
data  any  error  made  in  measuring  the  echo  height  chosen  as  the  basis  for 
comparison.   However,  Table  I  shows  that  essentially  the  same  results 
are  obtained  using  either  method.   The  difference  in  the  methods  is 
apparent  in  the  different  values  of  residual  attenuation  obtained.   Taking 
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the  value  at  the  lowest  temperature  reached  for  comparing  echo  heights 
required  that  the  residual  attenuation  be  determined  from  this  value 
rather  than  from  an  average  of  low  temperature  values  as  was  done 
using  semi- log  plots. 
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TABLE  1 
RESULTS  OF  RESATT  COMPUTATIONS 


Run 


13  Nov  64 


5  Jan  65 


5  Jan  65 


Attenuation  from  Semi- log  Plots 

Frequency       Normal      Residual    2£(0) 
Mcs/sec    Attenuation   Attenuation    KT 


30 


30 


10 


OL   db/cm 
n 

a   , db/cm 

<_ 

0.310 

3.46 

3.092 

0.320 

3.58 

0,330 

3.62 

0.525 

2.96 

3.350 

0.565 

3.26 

0,605 

3.49 

0,660 

2.73 

1.376 

0,700 

3.25 

0.740 

3.63 

Results 


CC  ,  db/cm 


0.320 


0.565 


0.700 


2€(0) 

KT 


3.58 


3.26 


3.25 


Attenuation  from  Comparison  of  Echo  Heights 


Run 

Frequency 
Mcs/sec 

Normal 

Attenuation 

CC  db/cm 
n 

Residual 

Attenuation 

OC    ,  db/cm 
o 

2€(0) 

KT 
c 

Result 

s 

a  , db/cm 
o 

2€(0) 

KT 
c 

0.290 

3.45 

13 

Nov 

64 

30 

3.092 

0.300 
0.310 

0.525 

3.58 
3.65 

2.75 

0.300 

3.58 

5 

Jan 

65 

30 

3.350 

0.565 
0,605 

0.700 

2.99 
3.57 

3.01 

0.585 

3.30 

5 

Jan 

65 

10 

1.376 

0.715 
0.730 

3.33 

0.715 

3.33 
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DISCUSSION  OF  RESULTS 

The  value  obtained  for  the  transition  temperature,  T  =  0.817  K, 

is  lower  than  those  of  other  investigators.   However,  those  other 

experiments  were  determinations  of  the  zero  magnetic  field  transition 

temperature  for  zinc  and  show  considerable  variation.   In  heat  capacity 

measurements,  Phillips  determined  T  =  0,825  K.   Geballe  and  Matthias 

r  c 

in  work  on  isotope  effect  measured  T  =  0,855  K  and  Cochran  and 

11  C 

Mapother  *  in  critical  field  experiments  found  T  ■  0,875  K,   If  a 

parabolic  relationship  between  reduced  magnetic  field  and  reduced 

temperature  is  assumed 

.2 

c 

where  H   is  the  critical  magnetic  field  above  which  zinc  will  not 
c 

become  superconducting,  then  the  transition  temperature  is  reduced 
about  five  millidegrees  for  H  =50  gauss  and  H  =  0,5  gauss.  Thus 
T  =  0,817  K  agrees  favorably  with  the  lowest  value  quoted  above. 

Table  II  gives  the  results  of  several  investigators  who  have 
measured  the  zero  degree  energy  gap  in  zinc.   There  have  been  no  ultra- 
sonic attenuation  measurements  along  the  c-axis  other  than  this  one. 
The  result  of  this  work  is  higher  than  the  value  found  by  Biondi  et  al 
in  a  microwave  absorption  experiment  in  which  there  were  "violent 

oscillations"  in  the  data  for  hv/kT   less  than  2,7,   The  measurement 

c 

by  Zemon  and  Boorse  can  only  be  taken  as  an  average  over  all  directions 
in  the  crystal,   Bohm  and  Horowitz  found  the  attenuation  so  great 
along  the  c-axis  that  they  were  unable  to  make  any  measurements.   How- 
ever, their  measurement  of  the  energy  gap  in  two  directions  in  the  base 
plane  and  the  measurements  of  the  energy  gap  along  the  c-axis  by 
Biondi  and  this  work  show  anisotropies  in  the  energy  gap.   These  results 
also  indicate  that  the  gap  along  the  c-axis  is  less  than  that  predicted 
by  the  BCS  theory  and  also  is  less  than  that  for  directions  in  the 
base  plane. 
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TABLE  II 
ZERO  DEGREE  ENERGY  GAP  IN  SUPERCONDUCTING  ZINC 


2  (0)/kT 


3.8  +  0.2 


3.4  +  0.2 


3.00  +  0.15 


3.01  +  0.09 


3.36  +  0.13 


Direction 

T 
c 

0.85 

Method 

Investigator 

[l2lo] 

Ultrasonic 
Attenuation 
233  Mcs/sec 

Bohm  and 
Horowitz 

[  ioTo] 

0.85 

Ultrasonic 
Attenuation 
233  Mcs/sec 

Bohm  and  . 
Horowitz 

[oooi] 

0.838 

Microwave 
Absorption 
15-100  Gcs/sec 

BiondJ 
et  al 

Polycrystalline 
Sample 

0.825 

Microwave 
Absorption 
50-75  Gcs/sec 

Zemon  and 
Boorse 

[oooi] 

0.817 

Ultrasonic 

Attenuation 

10  6c  30  Mcs/sec 

Present 
Work 
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A  more  precise  determination  could  be  made  if  lower  temperatures 
could  be  reached.   The  lowest  temperature  reached  in  this  work  was 
0.320  K;  i.e.,  0.38  t.   Thus  to  determine  the  energy  gap,  the  experi- 
mental data  had  to  be  extrapolated  some  40%  on  the  reduced  temperature 
scale . 

The  graphs  of  reduced  attenuation  vs  reduced  temperature  can  be 
fitted  with  a  BCS  curve  quite  well  for  low  t.   However,  for  t  approaching 
1.00,  the  experimental  points  lie  above  the  BCS  curve.   Bohm  and  Horowitz 
don't  report  the  comparison  of  their  points  with  the  BCS  curve.   How- 
ever, investigators  of  other  elements  have  found  just  the  opposite; 
i.e.,  points  for  higher   t   lie  below  the  BCS  curve.   For  Niobium, 
Weber   found  the  effect  to  be  frequency  dependent.  With  T  =  9,15  K, 
the  BCS  curve  for  the  measured  value  of  the  energy  gap,  3,63  kT   , 
fit  the  experimental  points  for  220  Mcs/sec  longitudinal  wave  along 
the  [100]  axis.   However,  for  30  Mcs/sec  the  points  lay  below  the  BCS 
curve.   He  also  found  an  isotropic  result  for  the  energy  gap  in  the 
three  crystallographic  directions  and  suggested  that  any  anisotropy  in 
the  gap  may  have  been  smeared  out  by  impurity  scattering.   For  thallium, 
T  =  2.38  K,  Saunders  and  Lawson   found  that  attenuation  of  longitudinal 
waves  decreased  more  sharply  below  T  than  predicted  by  the  BCS  theory 

and  that  the  anisotropy  in  the  energy  gap  ranged  from  (3.76  -  4.10)kT  . 

c 

These  and  similar  results  on  other  elements  are  in  contrast  to  the 
findings  of  this  work  which  should  be  investigated  further.   One  basic 
assumption  tacit  in  all  these  computations  was  that  the  residual  atten- 
uation was  temperature  independent.   Though  measurement  above  the 
transition  temperature  bore  this  out,  the  attenuation  in  the  normal 
state  for  temperatures  less  than  the  zero  magnetic  field  T   should  be 
measured.   This  might  be  done  by  placing  the  sample  in  a  magnetic 
field  to  destroy  the  superconducting  state. 

Finally,  measurements  at  higher  frequencies  should  be  carried  out. 

Bobetic   has  carried  out  the  calculations  giving  the  ration  of  Oi   /a 

s   n 
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according  to  the  BCS  theory  for  an  isotropic  model  at  arbitrary 
frequencies  and  temperatures  and  predicted  a  discontinuity  in  attenu- 
ation at  frequencies  corresponding  to  the  energy  gap,   This  would 

2 

require  frequencies  of  the  order  10  Gcs/sec, 

In  this  work  attempts  were  made  to  make  measurements  at  50  and 
70  Mcs/sec,   However,  so  much  power  was  needed  to  get  an  echo,  that 
the  sample  heated  excessively.   Using  a  shorter  sample  would  help 
solve  this  problem . 
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APPENDIX  I 
DATA  ANALYSIS 

Analysis  begins  by  measuring  the  echo  heights  on  each  picture. 
This  is  most  accurately  done  on  an  X-ray  film  reader  to  +  0.01  mm.   These 
heights  are  then  plotted  on  semi-log  graph  paper  vs  echo  number;  one 
graph  per  picture.   The  slope  of  this  graph  is  absolute  attenuation  per 
echo. 

The  basic  problem  in  determining  the  absolute  attenuation  in  the 
temperature  range  attained  (4,2  K  to  approximately  0.320  K)  is  that  the 
attenuation  is  so  great  from  slightly  below  the  superconducting  transi- 
tion temperature,  Ts   to  4<,2  K  that  one  usually  has  only  one  echo. 

c 

Below  T  ,  the  attenuation  drops  rapidly  and  enough  echoes  appear  to 
obtain  the  attenuation  graphically. 

A  typical  plot  of  absolute  attenuation,  GC,    vs  temperature,  T,  is 
as  follows: 


T 
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where 

a  =  absolute  attenuation 

CU  -  residual  attenuation 

o 

d     -   attenuation  due  to  conduction  electrons  in  the  superconducting 
state 

Q:  =  attenuation  due  to  conduction  electrons  in  the  normal  state 
N 

a.    -   attenuation  determined  from  raw  data 

1 

T  -  Temperature  -   K 
T  =  Superconducting  transition  temperature 
ACC  =  Change  in  attenuation  between  the  superconducting  and  normal  states 

The  slope  of  the  semi- log  plots  are  in  echo   .   The  attenuation  in 
db/cm  is  given  as  follows: 

/ju/   s    Slope  of  Semi-log  Plot  x  8.6858 

Attenuation  (db/cm)  =  r — c : : -=— ; 

v     '    2  x  sample  length  in  cm 

We  assume  an  exponential  decay  of  echo  heights . 

A  =  A'e 

Let 

t-Vi 

A.  -  amplitude  of  n   echo  in  the  superconducting  state 

A  =  amplitude  of  n   echo  in  the  normal  state 
N 

A'  =  amplitude  of  the  initial  pulse 

n  -  echo  number 

CUT  =  attenuation  in  the  normal  state  from  T   to  4.2  K;  a  constant 
c  c 

as  verified  experimentally 

OL.    =   temperature  dependent  attenuation  in  the  superconducting  state 

L  =  sample  length  in  cm 

Then  rv   •  n 

A  =  A'e    ° 
AN  -  A  e 

-a.  -n 

A.  =  A'e 

l 

Ln  A„  =  At   A8  -  n« 
N  T 

c 

Ln  A.  -  /n  A1  -  n  a. 
i  l 

A. 

in(^)-  n  (a  -  ai)  =  -nAa 

N        c  -31- 


,...     N   8.6858  D     ,Ai  . 
or    aT  (db/cm)  =  j-^ Xn   (— )  +  QL 

c  N 

To  determine  OC       from  the  single  echo  pictures  above  T  ,  one  com- 
*-c  c 

putes  OC     ,  with  n  =  15  for  as  many  point  by  the  above  equation  as  one 

Tc 
can  from  the  semi-log  plots.   The  best  horizontal  straight  line  fit  to 

a  plot  of  these  points  gives  a  .   Then  using  the  equation 

Tc 


,..,     s        8.6858  /  .   TI  v 


i     "« 

gives  the  attenuation  for  the  remaining  points  below  T  where  there 

are  not  enough  echoes  to  determine  the  attenuation  from  the  semi- log 

plots. 

Extrapolating  the  plot  of  a.  and  a       vs  T  to  T  -  0  will  give  an 

1      Tc 

estimate  of  OC    .      Where  the  curved  part  intersects  the  value  for  QL, 
o  Tc 

determines  T  . 

c  *  as 

The  quantity  desired  is  the  reduced  attenuation  Oi     =  ,  the 

% 
ratio  of  attenuation  due  to  conduction  electrons  alone  in  the  super- 
conducting and  normal  states  vs  t  =  T/T  ,  the  reduced  temperature. 
Subtracting  from  a.  and  _   the  quantity  OC     which  is  attributed  to 
impurity  scattering.,  end  lossess  transducer  losses,  etc.,  gives 

a,    a.   -  a 
*    S     i    o 
a  - 


aN    aTc  "  ao 


Subsequent  calculations  are  very  sensitive  to  the  value  Ct   ,  which 
is  not  easily  determined  since  one  must  extrapolate  approximately  40% 
on  the  reduced  temperature  scale  to  t  =  0  for  zinc. 

The  BCS  theory  gives 

*       2 

a 


^€(T)/kT+  l 

for  the  variation  of  a     for  the  longitudinal  attenuation  with  T,  where 

€  (T)  is  the  temperature  dependent  energy  gap  between  normal  and  super- 
conducting electrons  and  k  is  the  Boltzmann  constant.  What  we  want  is 

-32- 


e (0) .   The  above  expression  can  be  inverted  to  give 

,  kT 


[l/a*  -   i] 

This  is  the  equation  of  a  straight  line  passing  through  the  origin 

* 
when  2/a  »  1.   Hence,  the  slope  at  t  =  0  for  a  plot  of  y    vs  t  gives 

the  reciprocal  of  the  zero  degree  superconducting  energy  gap  in  units 

of  kT  . 
c 

Data  can  also  be  analyzed  by  comparing  the  echo  height  of  the 
same  echo  on  each  of  the  pictures.   Then  for  a  given  n 

w  ...  ,   .    8,6858  a     ,   Ano  s 

a.(db/cm)  =  — Jfn    ( )  +  a   . 

i  2nL         A  .      o 

m 

One  must  have  some  reference  attenuation  here  denoted  OL     to  signify 

o 

the  attenuation  at  the  lowest  temperature  reached.  The  attenuation  can 
be  most  accurately  determined  at  this  point  since  the  maximum  number  of 
echoes  is  present. 
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APPENDIX   II 
PROGRAM  RESATT 

The  raw  data;  i.e.a  pictures,  must  be  first  analyzed  in  the  manner 
outlined  in  Appendix  I.   Then  the  computer  can  be  utilized  to  compute 
the  reduced  attenuation  and  the  function  Y,  used  to  determine  the  zero 
degree  superconducting  energy  gap,  and  plot  these  values  separately  as 
a  function  of  the  reduced  temperature.   Program  RESATT,  written  in 
Fortran  60  language  and  compatible  with  the  1604  CDC  computer,  will 
accomplish  these  calculation.   The  primary  purpose  of  this  program  is 
to  determine  the  residual  attenuation  and  to  find  the  zero  degree  energy 
gap.   The  residual  attenuation  will  be  that  value  which  when  subtracted 
from  the  experimental  values  of  attenuation  will  make  Y,  for  t  less  than 
0.4,  plot  as  a  straight  line  intercepting  the  origin. 

For  each  set  of  data  the  program  provides  for  the  point  plotting 
of  computations  based  on  three  chosen  values  of  the  reduced  attenuation, 
each  set  plotted  on  a  separate  graph.   On  each  of  these  plots  a  theo- 
retical curve  based  on  the  BCS  expression,  using  a  chosen  zero  degree 
energy  gap>  can  be  graphed.   This  is  computed  from  values  of  £(T)/6(o) 

vs  t,  the  ratio  of  the  temperature  dependent  energy  gap  to  the  zero 

18 
degree  energy  gap  versus  reduced  temperature  as  follows   : 

*        2 
a  = 


C(I)EGAP 
2t 

X  +  1 

Presently  44  values  of  C(I)  are  used.   More  can  be  computed  from  the 
reference. 

The  titles  which  must  be  changed  for  each  set  of  data,  identifies 
the  experimenter ,  date  on  which  datawere  taken,,  frequency,  and  computer 
project  number.   These  are  within  the  program  itself.   The  program  notes 
follow. 

Code 

TTR(I)  -  Reduced  temperature,  t,  for  theoretical  BCS  plot 

•k 

ATS  (I)  -  Reduced  attenuation,,  C£  ,  for  theoretical  BCS  plot 

-34- 


•k  -1 

TY(I)  -  (ln(2/a  -1))   ,  for  theoretical  BCS  plot 

M  -  Number  of  data  points  for  theoretical  BCS  plot,  maximum  number,  100 

EGAP  -  Zero  degree  energy  gap,  2  £ (0) ,  (no  kT  ) 

C(I)  -  6(TTR(I)/6(0) 

N  -  Number  of  experimental  values  of  T  &  (X.    (up  to  60) 

TC  -  Experimental  value  of  Transition  Temperature 

AN  -  Experimental  value  of  a  =  QL, 

n    Tc 

AZ(0)  -  Three  experimental  values  of  OL 

T(I)  -  N  experimental  values  of  temperature 

A(I)  -  N  experimental  values  of  attenuation 

RT(I)  -  N  values  of  reduced  temperature 


AS(1,I) 
AS(2,I) 
AS  (3, 1) 


N  values  of 

*   a  -  af 

a  =  — 


a 


n 


a 


( AZ(0)  -  small 
AZ(0)  -  medium 
AZ(0)  -  large 


AY(1,I) 
AY(2,I) 
AY(3,I) 


N  values  of 

1 


y  = 


U 


a. 


-1 


( AZ(0)  -  small 
AZ(0)  -  medium 
AZ(0)  -  large 


Plots:  Two  sets  of  point  plots 

First:  Three  plots  of  a     vs  t, 


Experimental  Points  AS(1,I)  vs  RT(I)  for  AZ(0)  small  -  □ 
Experimental  Points  AS(2,I)  vs  RT(I)  for  AZ(0)  medium  -  O 
Experimental  Points  AS(3SI)  vs  RT(I)  for  AX(0)  large  -  2\ 

Each  will  have  the  theoretical  BCS  curve  ATS(I)  vs  TTR(I)  plotted 
on  the  same  graph. 

Second:  Three  plots  of  Y  vs  t, 

Experimental  Points  AY(1,I)  vs  RT(I)  for  AZ(0)  small  -  O 
Experimental  Points  AY(2,I)  vs  RT(I)  for  AZ(0)  medium  -  ^ 
Experimental  Points  AY(3,I)  vs  RT(I)  for  AZ(0)  large  -  A 

Each  of  these  will  have  the  theoretical  BCS  curve  TY(I)  vs  TTR(I) 
plotted  on  the  same  graph. 
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For  each  run  a  total  of  3  +  M  +  N  cards  are  needed  in  the  following 
order.   The  last  digit  of  each  imput  in  each  column  in  the  stated  order 
is  to  be  in  the  10,  20,  30,  etc.,  column  of  the  data  card. 

Number  of  cards  Imput 

1  M,  EGAP 

1  -  M  Card  No.,  TTR(I) ,  C(I) 

1  N,  AN,  TC 

1  -  N  Card  No.,  T(I),  A(I) 

1  AZ(1),  AZ(2),  AZ(3) 

Following  is  a  complete  print  out  of  RESATT  with  the  card  deck 
composed  for  the  data  of  13  Nov  64,   Card  000-197  comprise  the  program 
and  cards  200-243  are  the  Mtlhlschlegel  values.   As  outlined,  the  cards 
to  be  changed  for  a  different  set  of  data  would  be  57-61,  199,  and 
244-265. 
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APPENDIX  III 
PROGRAM  VARIGAP 

Once  the  best  value  of  the  residual  attenuation  and  a  determination 
of  the  zero  degree  energy  gap  has  been  made  from  RESATT,  VARIGAP  will 
compute  and  plot  the  theoretical  BCS  curves  of  reduced  attenuation  vs. 
reduced  temperature  and  the  function  Y  vs  reduced  temperature  for  three 
different  zero  degree  energy  gaps.   On  each  of  these  will  be  the  cor- 
responding experimental  points  for  one  set  of  data  and  the  chosen  re- 
sidual attenuation.   To  the  extent  that  they  differ  from  RESATT,  the 
program  notes  for  VARIGAP  follow. 

CODE 

EGAP(I)  -  Three  values  of  the  zero  degree  energy  gap  (no  kTc) 

AS      -  One  value  of  the  residual  attenuation 

ATSCljI)^    M  theoretical  values  of  (EGAP  -  small 

ATS (2, I) 7    reduced  attenuation  from      for   S EGAP  -  medium 
ATS(3,I)J     the  BCS  expression  (eGAP  -  large 

TY(1,I)]     M  theoretical  values  of  ^EGAP  -  small 

TY(2,I)>     Y  from  the  BCS  expression     for   <  EGAP  -  medium 
TY(3,I)J  (_EGAP  -  large 

AS (J)    -  N  experimental  values  of  reduced  attenuation 
AY(J)    -  N  experimental  values  of  Y 

PLOTS :   Two  sets  of  curve  plots 

First:  Theoretical  reduced  attenuation  vs  reduced  temperature 

ATS (1,1)  vs  TTR(I)  for  EGAP  -  small  -  □ 
ATS(2,I)  vs  TTR(I)  for  EGAP  -  medium  -<Q 
ATS (3, I)  vs  TTR(I)  for  EGAP  -  large  -  A 

Each  of  these  will  have  point  plotted  on  them  the  N  experimental 
values  of  AS(J)  vs  T(J) . 

Second:  Theoretical  Function  Y  vs  Reduced  Temperature 

TY(1,I)  vs  TTR(I)  for  EGAP  small 
TY(2,I)  vs  TTR(I)  for  EGAP  medium 
TY(3,I)  vs  TTR(I)  for  EGAP  large 
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Each  of  these  will  have  point  plotted  on  them  the  N  experimental 
values  AY(J)  vs  T(J) . 

For  each  run  a  total  of  3  +  M  +  N  cards  are  needed  plus  the  correct 
title  cards .   The  last  digit  of  each  imput  in  each  column  in  the  stated 
order  is  to  be  in  the  10,  20,  30,  etc.,  column  of  the  data  card. 

Number  of  Cards  Imput 

1  M,  EGAP(low),  EGAP(medium),EGAP( large) 

1  -  M  Card  No.,  TTR(I),  C(I) 

1  N,  AN,  TC 

1  -  N  Card  No.,  T(J),  A(J) 

1  AZ 

Following  is  a  complete  print  out  of  VARIGAP  with  the  card  deck 
composed  for  the  data  of  13  Nov  64 .   Cards  000  -  188  comprise  the  program, 
As  outlined,  the  cards  to  be  changed  for  a  different  set  of  data  would 
be  57  -  61,  190,  and  235  -  257. 
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APPENDIX  IV 
PROGRAM  ZEROGAP 
This  program  provides  for  a  least  square  fit  of  experimental  points 
to  the  BCS  function 

s  2 


by  minimizing  the  square  of  the  vertical  distances  from  the  curve  to 
the  experimental  points.   This  is  done  using  two  parameters;  the  re- 
sidual attenuation  and  the  zero  degree  energy  gap.   The  output  of  the 
program  is  the  value  of  these  two  parameters,,   Since  in  the  BCS  ex- 
pression a  and  a  involve  only  the  attenuation  due  to  electron-phonon 
s      N 

interaction,  all  other  causes  of  attenuation  must  be  subtracted  out. 
Thus, 


a1 -a 
o 

2 

N      o 

„  kT       i 

or 


a'  =  a  +  (a„-a  ) 


o   N  N  o'  £(T) 

where  a'  and  QL.  would  be  on  a  BCS  curve  shifted  an  amount  a     due  to 
N  o 

residual  attenuation. 

What  is  desired  is  to  minimize  the  residue 

N 

r  =  Y.      (a-a')2 
i=l 

where  a  is  an  experimental  measurement  of  attenuation  and  N  is  the  number 

of  experimental  points. 

Let  the  following  quantities  be  defined  which  will  be  used  in  the 

program. 

CH  -  A(I)  -  Experimental  measurement  of  attenuation 

t  -  RT(I)-  The  reduced  temperature  corresponding  to  the  temperature  T(I) 
at  which  A(I)  was  measured. 
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OL   -   AZ  -  Residual  attenuation  -  one  of  two  parameters  to  be  determined 
o 

C(I)-US(I)  -  The  Miihlschlegel  values  of  £(T)/  6(0)  vs  t. 


By  use  of  an  interpolating  program,  the  values  of  C(I)  were  obtained 
from  the  Miihlschlegel  values  for  0.191  C   t   <.  1.000  in  increments  of 
0.001  and  incorporated  in  the  program.   In  the  printout  of  the  program 
following  they  are  cards  200  -  910,  where  for  example,  200  means 
t=0.200  and  the  number  in  the  10th  column  is  C(I)  for  the  t. 
C(I)/t  -  UC(I) 

The  equation  for  CC*    can  then  be  written 


"'  =  ao  +   <Vao>    c(I)-  2g(0) 


£ 


2t 


+  1 


Now  let  FER(I) 


C(I)-  2  €(0) 

2t  EXP 

Z  +i 


UC(I) -EGAP 
2 


+  1 


where  EGAP  is  the  zero  degree  energy  gap  =  2  6(0) 
to  be  determined.   Then  the  residue  becomes 
N 


the  other  parameter 


I 

i=l 


A(I)-AZ-(AN-AZ)2FER(I) 


From  this  an  explicit  expression  for  AZ  can  be  gotten  from  the 
condition 


«2R 

5AZ 


-  0, 


Carrying  out  the  steps  gives 
N 


AZ 


£     [(2ANFER(I)-A(I)1    2FER(I)-l] 
£x   I  2FER(I)-l]  [2FER(I)-l] 


in  which  terms  cannot  be  canceled  since  they  are  under  the  summation 
signs.   In  this  expression  only  FER  depends  on  EGAP. 
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5R 

The  other  condition  for  minimizing  the  residue  is   ^      =  0. 

e»  EGAR 

Carrying  out  the  steps  gives 

N 
gP~£p  -  2  Z   [A(I)-AZ-2(AN-AZ)FER(I)] 

X   [2(AZ-AN)  -~  FER(FER-l)] 

So  what  is  needed  is  to  minimize   £  f(2ANFER(I) -A(I)) -AZ(2FER-1)1  UC 

X  [  FER(FER-l)!  I"aN-Az1 

subject  to  the  above  value  of  AZ»   The  remainder  of  the  program  notes 
are  the  following: 

ONE (I)  -  2FER(I)  -  1 
TWO(I)  =  2ANFER(I)  -  A(I) 

The  imputs  needed  are  the  best  guess  of  EGAP;  AN,  the  experimental 
determination  of  the  attenuation  in  the  normal  state;  TC,  the  transition 
temperature;  and  the  N  data  points  of  A(I)  vs  T(I) . 

One  card  contains  N,  AN,  TC,  and  EGAP,  the  last  digit  of  each  is  in 
column  10,  20,  30  and  40  respectively.   Next  follow  the  N  data  cards  of 
A(I)  vs  T(I) ;  one  card  for  each  set  with  the  last  digit  of  the  card  number 
in  column  10,  that  of  T(I)  in  20,  and  that  of  A(I)  in  30.   A  print  out  of 
the  program  follows  at  the  end  of  the  appendix .   The  imput  cards  follow 
the  blank  card  which  must  be  added  to  the  end  of  the  program  cards. 

Table  III  gives  the  results  of  program  Zerogap  and  compares  them 
to  those  of  program  Resatt.   Figs.  15-26  give  the  plots  of  the  BCS  curve, 
using  the  values  of  residual  attenuation  and  the  energy  gap  determined 
from  Zerogap  and  the  point  plots  of  the  corresponding  experimental  data. 

It  is  evident  that  the  program  does  what  it  is  designed  to  do. 
However,  the  values  of  residual  attenuation  are  markedly  different  in  the 
30  mcs/sec  data  from  that  obtained  by  extrapolation  of  the  curve  of 
attenuation  vs  temperature.   For  the  30  mcs/sec,  5  Jan  65,  semi-log  data 
the  value  obtained  is  unrealistic  in  that  it  is  greater  than  most  of  the 
data  points o   Only  the  attainment  of  lower  values  of  reduced  temperature, 
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t,  can  resolve  the  question  as  to  whether  or  not  the  value  of  residual 
attenuation  obtained  from  Zerogap  is  realistic . 

The  program  seems  to  work  best  when  there  are  data  points  over 
the  whole  range  of  t;  e.g.,  30  Nov  64  and  5  Jan  65,  10  mcs/sec  data. 
Also,  when  the  points  are  not  too  scattered  as  in  the  30  Nov  64  data, 
the  results  are  reasonable.   The  values  of  (2„85  and  2.82)kT  obtained 
for  the  energy  gap  for  these  data  are  reasonable  bearing  in  mind  the 
previous  observation  that  for  higher  t  the  experimental  points  lie  above 
the  3.50  kT  BCS  curve.   As  one  can  see  the  lower  the  value  of  the  energy 
gap,  the  higher  the  BCS  curve. 

This  program  can  be  of  use  in  interpreting  data  if  one  realizes 
that  the  analysis  depends  completely  on  the  BCS  theory  and  the  simplifying 
assumptions  used  in  obtaining  the  expression  for  the  reduced  attenuation 
and  masks  the  possibility  of  other  phenomenon  taking  place. 
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TABLE  III 
COMPARISON  OF  ZEROGAP  AND  RESATT  RESULTS 


RUN 


FREQUENCY 
Mcs/sec 


TYPE 


ZEROGAP 


RESIDUAL   ENERGY 
ATTENUATION    GAP 
db/cm     2€(0)/kT 


RESATT 


RESIDUAL   ENERGY 
ATTENUATION   GAP 
db/cm    2€(0)/kT 


13 

Nov  64 

30 

Semi-log 

0.1772 

2.8544 

0.320 

3.58 

13 

Nov  64 

30 

Echo  Heights 

0.1507 

2,8171 

0.300 

3.58 

5 

Jan  65 

30 

Semi- log 

0.8755 

7.2968 

0.565 

3.26 

5 

Jan  65 

30 

Echo  Heights 

0,4387 

2 . 5304 

0,585 

3.30 

5 

Jan  65 

10 

Semi-log 

0.7166 

3.5941 

0.700 

3.25 

5 

Jan  65 

10 

Echo  Heights 

0.6097 

2.1181 

0.715 

3.33 
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APPENDIX  V 
PROGRAM  AVERGAP 

This  program  provides  for  a  least  squares  fit  to  the  BCS  function 
using  one  parameter;  the  zero  degree  energy  gap  with  the  residual  atten- 
uation specified.   Thus  the  program  and  the  program  notes  are  the  same  as 
those  of  Program  Zero-gap,  except  that  the  residual  attenuation  is  part 
of  the  imput  data  rather  than  being  derived  from  an  analytic  expression. 
The  program  will  compute  the  zero  degree  energy  gap  which  minimizes  the 
square  of  the  vertical  distance  from  the  experimental  points  to  the  BCS 
curve  for  three  different  values  of  residual  attenuation,   Thus,,  the 
expression 

V     |"(2ANFER(I)  -  A(I))  -  AZ(2FER(I)  -1)   UC  [  FER(FER-ll  [aN-Az] 
i=l 

is  minimized  for  each  given  AZ . 

The  type  of  data  cards,  their  format,  and  order  of  imput  are  the  same 
as  for  Program  Zerogap,  with  the  exception  that  one  card  need  be  added. 
This  is  the  last  card  and  contains  the  three  values  of  residual  atten- 
uation; the  last  digit  of  each  to  end  in  column  10,  20,  and  30  respectively. 
These  values  are  to  be  determined  as  in  Program  Resatt;  i.e.,  from  an  ex- 
trapolation to  T  =  0  on  a  plot  of  the  experimental  data  of  attenuation  vs 
temperature. 

Table  IV  gives  the  results  of  Program  Avergap  and  also  a  summary  of 
all  the  computer  calculations.   It  is  seen  that  the  results  of  Avergap 
are  more  consistent  with  those  of  Resatt  than  are  those  of  Zerogap  and 
again  reflect  the  fact  that  the  experimental  points  for  higher  t  lie  above 
the  BCS  curves  used  in  Resatt.   The  results  for  the  data  of  5  Jan  65,  30 
Mcs/sec  still  seem  unreliable  with  such  a  high  energy  gap.   This  is  due  to 
the  fact  that  there  are  data  points  for  low  t  only. 
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